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ABSTRACT 

A  highly-scalable,  Bayesian  approach  to  the  problem  of  performing  multi-source  data  fusion  and  target 
tracking  in  decentralized  sensor  networks  is  presented.  Previous  applications  of  decentralized  data  fusion 
have  generally  been  restricted  to  uni-modal/uni-source  sensor  networks  using  Gaussian  based  approaches, 
such  as  the  Kalman  or  information  fdter.  However,  with  recent  interest  to  employ  complex,  multi¬ 
modal/multi-source  sensors  which  potentially  exhibit  observation  and/or  process  non-linearities  along  with 
non-Gaussian  distributions,  the  need  to  develop  a  more  generalized  and  scalable  method  of  decentralized 
data  fusion  is  required.  The  probabilistic  approach  featured  in  this  work  provides  the  ability  to  seamlessly 
integrate  and  efficiently  fuse  multi-source  sensor  data  in  the  absence  of  any  linearity  and/or  normality 
constraints.  The  proposed  architecture  is  fully  decentralized  and  provides  a  methodology >  that  scales 
extremely  well  to  any  growth  in  the  number  of  targets  or  region  of  coverage.  This  paper  will  illustrate  that 
our  multi-source  data  fusion  architecture  is  capable  of  providing  high-precision  tracking  performance  in 
complex,  non-linear/non-Gaussian  operating  environments.  In  addition,  we  will  also  show  that  our 
architecture  provides  an  unprecedented  scaling  capability  for  decentralized  sensor  networks  as  compared  to 
similar  architectures  which  communicate  information  using  particle  data,  Gaussian  mixture  models  or  Parzen 
density  estimators. 


1.0  INTRODUCTION 

In  the  field  of  multi-sensor  data  fusion,  decentralized  data  fusion  has  become  an  attractive  alternative  to 
centralized  data  fusion  primarily  due  to  the  inherent  robustness  and  scalability  features  that  decentralized 
architectures  offer.  In  its  most  primitive  form,  a  decentralized  sensor  network  involves  processing  capability 
at  each  sensor  -  eliminating  the  need  and  subsequent  vulnerability  of  a  central  processing  node  -  along  with 
the  capacity  for  each  sensor  to  efficiently  communicate  its  information  to  neighboring  sensors  without 
requiring  any  knowledge  of  the  network  topology  [1]. 

To  date,  the  majority  of  fielded  implementations  utilizing  decentralized  data  fusion  have  relied  on 
linear/Gaussian  assumptions  and  the  Kalman/information  filter  [2, 3, 4, 5, 6].  Even  though  such  systems  have 
produced  impressive  results,  the  natural  desire  to  utilize  a  wider  mixture  of  more  complex  sensor  types  - 
potentially  exhibiting  observation  and/or  process  non-linearities  along  with  non-Gaussian  distributions  -  has 
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generated  a  need  for  more  generalized  information  fusion  techniques.  A  variety  of  methods  have  been  applied 
to  the  problem  of  non-linear/non-Gaussian  decentralized  data  fusion  in  which  the  majority  of  such  methods 
have  been  focused  on  particle  filters  [7,8],  Gaussian  mixture  models  or  Parzen  density  estimators  [9,10]  or 
some  combination  of  the  two  [11,12,13].  All  of  the  aforementioned  techniques  are  capable  of  performing 
generalized  estimation,  however,  not  all  techniques  lend  themselves  well  to  decentralized  applications  as  a 
result  of  the  scalability  concerns  they  inherently  generate. 

This  paper  presents  a  highly-scalable,  Bayesian  approach  to  the  problem  of  performing  generalized,  multi¬ 
source  data  fusion  and  target  tracking  in  decentralized  sensor  networks.  A  brief  background  of  linear 
decentralized  data  fusion  is  first  provided  in  Section  2,  and  performance  results  related  to  currently  fielded 
linear  decentralized  data  fusion  implementations  are  provided  in  Section  31.  Section  4  introduces  a  proposed 
approach  to  performing  probabilistic  decentralized  data  fusion,  and  corresponding  simulated  performance 
results  are  presented  in  Section  52.  Section  6  summarizes  the  paper  and  provides  direction  for  future  efforts. 


2.0  LINEAR  DECENTRALIZED  DATA  FUSION 

2.1  Data  fusion  in  sensor  networks 

When  formulated  from  a  Bayesian  standpoint,  the  data  fusion  problem  is  to  recursively  compute  the  posterior 
probability  density  function  (PDF)  of  the  state  of  interest  jct  at  time  step  k ,  using  the  measurement  set  Zk  =  {zj, 
j=l,2,...,k };  this  can  be  calculated  using  Bayes  rule 

p(Xtlz*)=ZliM£ZiZZ2  (2.U) 

P(?k\Zk~l) 

where  p(zk\xk)  is  the  measurement  likelihood  and  p(xk\Zk^)  is  the  predicted  PDF  at  time  k  computed  according 
to  (2.1.2).  The  denominator  of  (2.1.1)  acts  as  a  normalisation  factor. 

p(xk\Z  )  =  |  p(xk\xk_i  )p(xk_x\Z  k~l  )dxk_x  (2. 1 .2) 

When  the  underlying  target  distribution  and  likelihood  are  assumed  to  be  Gaussian,  equations  (2.1.1),  (2.1.2) 
reduce  to  the  Kalman  filter  [14]  and  its  Information  filter  equivalent. 

To  allow  the  global  state  estimate  to  be  computed  in  a  decentralised  system  -  where  only  communication  with 
nearest  neighbours  is  allowed  -  nodes  must  communicate  their  own  data  and  data  received  from  other  sensors. 
Care  must  be  taken  to  prevent  nodes  receiving  data  along  multiple  paths.  If  not  correctly  handled  the  re¬ 
fusion  of  information  in  this  way  gives  rise  to  the  rumour  propagation  problem  which  can  seriously  degrade 
the  quality  of  fusion  [15].  It  is  also  desirable  that  each  node  should  not  have  to  record  the  provenance  of 
every  piece  of  information  generated  by  any  sensor.  For  tree  connected  communication  networks  the  solution 
to  this  problem  is  the  Channel  filter  [16].  The  Channel  filter  is  a  probability  distribution  conditioned  on  the 
common  information  held  between  a  pair  of  nodes.  When  Channel  filters  are  used  to  identify  redundant 
information,  only  posterior  distributions  need  to  be  transmitted  between  nodes.  This  allows  the  computation 
of  the  global  posterior  without  sending  measurements  or  provenance  information. 


1  Sections  2  and  3  are  provided  by  BAE  Systems  -  Advanced  Technology  Centre. 

2  Sections  4  and  5  are  provided  by  BAE  Systems  -  Advanced  Systems  &  Technology. 
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If  nodes  a  and  b  communicate  their  information,  the  updated  estimate  can  be  calculated  as  the  product  of  their 
distributions  divided  by  the  common  information  [12]. 


p(Xk\ZakvjZbk)  = 


P(xk\Zl)p(xk\Zbk) 

p(xk\ZaknZbk) 


(2.1.3) 


When  the  underlying  distributions  are  Gaussian,  the  channel  filter  is  also  an  Information  filter.  Importantly, 
since  information  is  additive,  redundant  information  can  be  removed  by  subtraction.  A  discrete  time  model  of 
a  linear  system xk  =  [xu  x2,...,xn],  XjGtR,  is  given  by  (2.1.4). 


xk  =  Fkxk- 1  +  Gkwk  (2.1.4) 

zk=Hkxk+vk  (2.1.5) 


where  Fk  is  the  state  transition  model;  Gk  projects  the  noise  vector  into  the  state,  and  wk  is  a  zero-mean,  white 
noise  vector  with  covariance  Qk.  An  observation  of  the  state  is  made  according  to  (2.1.5)  where  vk  is  a  zero 
mean  white  noise  vector  with  covariance  Rk. 

For  such  a  system  the  Kalman  filter  produces  an  optimal  estimate  x{k  \  k )  of  the  true  state  x(k)  and  an 
associated  covariance  P(k  |  k)  =  t[(  x(k  \  k )  -  x(k))'(  x(k  \  k)  -  x(k))\  using  the  measurement  set  Zk.  The  Kalman 
filter  equations  can  be  written  in  information  form  by  introducing  the  identities  (2.1.6)  for  the  information 
vector  and  information  matrix  respectively. 


y(k\k)  =  P(k\kylk(k\k) 
Y(k\k)  =  P{k\kyl 


(2.1.6) 


The  Kalman  filter  update  in  information  space  is  a  simple  addition  of  the  predicted  information  matrices  and 
vectors  and  the  measurement  written  in  information  form  (2.1.7),  (2.1.8).  The  predicted  densities  can  be 
calculated  by  converting  the  predicted  state  space  (2.1.9)  estimate  to  information  space. 


y(k\k )  =  y(k\k  - 1)  +  i(k) 
Y(k\k)  =  Y(k\k-\)  +  I(k) 

i(k)  =  HkRklzk 
I(k)  =  H'kR-lHk 


x(k  \k-\)  =  Fkx(k-\  |  k-\) 

P(k\k- 1)  =  FkP(k  - 1 1  k-  \)F'k  +  GkQkG'k 


(2.1.7) 


(2.1.8) 


(2.1.9) 


The  track-to-track  update  can  now  be  written  as  the  sum  of  information  vectors  and  matrices,  minus  the 
channel  filter  (2.1.10). 
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y{k  |  k)+  =ya{k\  k)  +  yb(k  \  k)-yab(k  \  k  -1) 

(2.1.10) 

Y(k\k)+  =Ya(k\k)  +  Yb(k\k)-Yab(k\k-\) 

The  channel  filter  is  then  updated  to  account  for  the  new  common  information  (2. 1 . 1 1). 

yab(k\k)= y(k\k)+ 

(2.1.11) 

Yab{k\k)  =  Y{k\k)+ 

Unfortunately,  the  channel  filter  can  only  be  used  in  single  connected  networks.  For  arbitrary  networks  the 
Covariance  Intersection  (Cl)  algorithm  can  be  used  [17].  The  Cl  algorithm  is  a  convex  combination  of  the 
information  vectors  and  matrices. 


y{k  |  k)+  =  coya  (k  \  k)  +  (1  -  co)yb  ( k  \  k) 

(2.1.12) 

Y(k\k)+  -coYa(k\  k)  +  (l-co)Yb(k\k) 

where  ms [0,1]  is  a  free  parameter  to  be  optimised,  typically  to  maximise  the  determinant  of  the  resultant 
information  matrix.  The  Cl  algorithm  gives  an  estimate  that  is  guaranteed  to  be  consistent  in  the  sense  that 
the  covariance  is  greater  than  the  expected  error.  Cl  thus  gains  additional  flexibility  at  the  expense  of 
optimality. 


3.0  LINEAR  DDF  -  FIELDED  PERFORMANCE  RESULTS 

3.1  DDF  Software  Development  Kit 

The  operations  of  a  Decentralised  Data  Fusion  (DDF)  node,  those  of  predict,  update,  communicate  and 
assimilate,  are  independent  of  the  underlying  representation  of  the  probability  distribution.  To  facilitate  the 
development  and  demonstration  of  novel  DDF  techniques,  BAE  Systems  ATC  has  exploited  this 
independence  and  developed  a  DDF  Software  Development  Kit  (SDK).  Illustrated  in  Figure  1,  the  DDF  SDK 
provides  an  overall  DDF  architecture  with  configuration,  sensor  and  communication  layer  interfaces  and 
support  functions  in  an  object-orientated  framework  written  in  C++.  Initially  the  SDK  utilised 
Kalman/Information  filters,  however  the  software  has  been  written  in  a  sufficiently  general  manner  that  any 
suitable  representation  such  as  particles  [12]  or  mixtures  of  Gaussians  [10]  can  be  added. 

3.2  Closed  loop  control 

The  DDF  Software  Development  Kit  (SDK)  has  been  used  by  BAE  Systems  ATC  to  develop  a  real  time 
closed  loop  sensor-to-actor  control  demonstration.  The  system  integrates  ATC  CCTV  cameras  (Figure  2a), 
Robovolc  mobile  robotic  capability  (Figure  2b)  and  Unattended  Ground  Sensors  (UGS)  (Figure  2c).  The 
scenarios  considered  were  used  to  demonstrate  the  advantages  of  DDF  within  a  heterogeneous  sensor  network 
in  an  end-to-end  system  that  encompassed  all  sections  of  the  OODA  loop  (Observe,  Orient,  Decide,  and  Act). 
The  actuator  was  the  Robovolc  robot,  the  decision  was  which  target(s)  to  intercept  and  when,  and  the  orient 
stage  involved  both  the  physical  orientation  of  the  robot,  and  the  sensor-to-target  allocation  for  steer-able 
sensors. 
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Figure  1:  DDF  SDK  is  a  data-fusion  software  architecture 


a 


b  c 

Figure  2:  CCTV  Camera  (a),  Robovolc  (b),  UGS  (c) 


The  scenarios  were  representative  of  a  number  of  military  scenarios  including  UAV  reconnaissance/ 
surveillance,  anti-submarine  warfare,  urban  asymmetrical  warfare  and  mine  counter  measures.  Figure  3a 
shows  the  experimental  set-up.  In  the  scenarios,  pedestrians  in  the  9K  car  park  give  rise  to  tracks  in  the 
Common  Operating  Picture  (COP).  Robovolc,  accessing  a  single  DDF  node  associated  with  a  CCTV  camera 
(marked  in  orange  in  Figure  3a),  defends  a  ‘sensitive’  region  of  the  car  park  against  intrusion  (e.g.  red  area  in 
Figure  3a).  Robovolc  defends  a  region  from  intrusion  by  physically  interposing  itself  in  the  path  of  targets 
that  appear  to  threaten  the  region.  Once  a  pedestrian  has  been  intercepted,  Robovolc  then  returns  to  guarding 
the  region.  Robovolc  moves  to  interception  points  supplied  to  it  by  software  that  sits  between  DDF  and 
Robovolc  (Figure  3b).  This  software  analyses  the  DDF  tracks  and  instructs  Robovolc  accordingly.  The 
instruction  to  move  is  triggered  by  a  target  crossing  an  outer  boundary  shown  in  green  in  Figure  3  a. 
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a  b 

Figure  3:  Operational  Setup  (a),  Software  Architecture  (b) 


The  UGS  sensor  (Figure  2c)  was  able  to  detect  targets  using  its  ultrasound  tripwire  while  running  the  SDK 
DDF  software  at  an  acceptable  rate  indicating  its  ability  to  participate  in  a  network  as  a  full  DDF  sensor  node. 

In  addition  to  the  DDF  tracking,  sensor  management  was  implemented  across  the  network.  Multi-sensor 
management  using  Decentralised  Decision  Making  (DDM)  algorithms  was  implemented  across  two  steer-able 
sensors  to  provide  a  globally  optimal  sensor  to  target  allocation.  The  manageable  sensor  was  mimicked  in 
software  by  only  processing  a  subset  of  the  cameras’  Field  Of  View  (FOV).  This  resulted  in  a  narrower  but 
steered  FOV.  The  sensor  management  software  controls  the  pointing  of  these  software-steered  cameras  with 
a  view  to  optimising  the  ability  of  Robovolc  to  defend  the  area.  The  results  are  drawn  from  analysis  of  the 
data  collected  during  two  experimentation  periods,  one  for  each  implementation.  The  measures  of 
effectiveness  are  drawn  from  the  SIAP  measures  [18].  The  relevant  SIAP  measure  is  completeness,  and  the 
relevant  system  metric  is  the  number  of  successful  (timely)  intercepts.  Figure  4  illustrates  two  still  images 
from  the  demonstration  video.  In  the  upper  right  and  left  corners,  the  steerable  camera  FOVs  are  shown  as 
light  grey  squares.  The  lower  left  comer  illustrates  a  plan  view  showing  the  tracked  target  means  and 
covariances. 


Figure  4:  The  UGS  was  used  as  a  trip  wire  to  cue  the  CCTV  FOV;  the  latter  are  the  light  grey  squares 
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4.0  PROBABILISTIC  DECENTRALIZED  DATA  FUSION 

When  applying  linear  decentralized  data  fusion  to  the  problem  of  target  tracking,  a  common 
observability/fusion  platform  in  position  space  is  typically  required  by  all  sensors  in  the  network  -  whether  it 
be  in  a  spherical  or  Cartesian  coordinate  system.  However,  when  multi-source  sensors  operating  in  disparate 
observability/fusion  spaces  are  used  to  populate  a  network,  a  common  platform  for  fusion  no  longer  exists, 
and  the  linear  decentralized  data  fusion  algorithms  must  be  abandoned  for  a  more  generalized  estimation 
approach. 

A  generalized,  probabilistic  approach  to  performing  multi-source  data  fusion  and  target  tracking  in 
decentralized  sensor  networks  is  illustrated  in  Figure  5.  The  proposed  architecture  provides  the  ability  to 
seamlessly  integrate  and  efficiently  fuse  multi-source  sensor  data  in  the  absence  of  any  linearity  and/or 
normality  constraints.  In  addition,  the  architecture  is  fully  decentralized  and  provides  a  methodology  that 
scales  extremely  well  to  any  growth  in  the  number  of  targets  or  region  of  coverage.  The  primary  components 
of  the  described  architecture  are  discussed  in  detail  in  the  subsequent  sections. 


* 


* 


Figure  5:  Probabilistic  Decentralized  Data  Fusion  Architecture 


4.1  Particle  Filters 

The  decentralized  probabilistic  approach  outlined  in  Figure  5  utilizes  particle  filters  to  estimate  the  local 
posterior  densities  for  each  sensor.  As  described  in  the  seminal  paper  by  Gordon  et  al.  [19]  along  with  others 
[14,20],  particle  filters  are  based  upon  point  mass  representations  of  probability  densities  which  can  be 
applied  to  any  state-space  model  in  the  absence  of  linearity  and/or  normality  constraints.  In  the  same  way  that 
a  Kalman  filter  provides  the  optimal  Bayesian  estimate  for  the  highly-restricted  linear,  Gaussian  environment, 
a  particle  filter  can  be  thought  of  as  approximating  the  optimal  Bayesian  recursion  without  having  to  impose 
any  linear  model  or  noise  distribution  constraints. 
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The  basic  idea  behind  the  particle  filter  is  that  the  posterior  density  of  interest  is  represented  by  a  set  of 
weighted  particles,  each  of  which  forms  an  independent  hypothesis  for  the  state  at  a  given  time.  If  the  weights 
are  chosen  correctly,  this  weighted  set  of  particles  becomes  representative  of  the  true  posterior  density  in  that 
expectations  of  the  true  posterior  can  be  made  arbitrarily  close  to  the  equivalent  expectations  of  the  set  of 
weighted  particles. 

In  order  to  describe  its  operation,  assume  that  a  particle  set  of  size  N  is  being  used  to  approximate  the 
posterior  density  of  a  state  vector,  jc*,  based  upon  the  set  of  all  available  measurements  Z*={zy,  j=l,2,...,k} .  At 
iteration  k~l,  the  posterior  density,  p(xk.I\Zk'1),  is  represented  by  a  set  of  particle  values  (x1  k.i,x2k.i,...,xNk.i)  and 
associated  weights  (w1  k.i,w2k.i,...,wNk.i),  At  the  next  iteration,  a  new  set  of  particle  values  (x1  k,x2k,...,x2'k)  and 
associated  weights  (w1  k,w2k,...,wNk)  which  characterize  the  posterior  density,  p(xk\Zk),  can  be  determined  as 
follows: 


xlrp(xk\xlkA) 

p(zk\xlk)wlM 
- _ 

i=l 

where  p(xk\x\.i)  is  the  transition  density  and  p(zk\xk)  is  the  likelihood  function.  Comparing  (4.1.2)  with  (2.1.1) 
it  can  be  seen  that  for  a  given  particle, 

<!=/>(*£  |ZM)  (4.1.3) 

wk=P(xk  lZ*)’  C4-1-4) 

As  a  result,  from  (2.1.1),  (4.1.3)  and  (4.1.4),  it  is  observed  that  the  weights  of  a  particle  filter  represent  the 
prior  and  posterior  densities  corresponding  to  each  particle’s  estimate  of  the  state  vector.  Consequently,  an 
overall  estimate  for  the  true  posterior  density  can  be  theoretically  obtained  from  the  particle  ensemble. 

4.2  Estimating  Local  Posteriors 

In  this  work,  range-only,  bearing-only  and  range/bearing  sensors  are  evaluated  for  multi-source,  decentralized 
data  fusion  in  which  particle  filters  are  used  to  estimate  the  local  posteriors  in  range/bearing  space  for  all 
sensor  types.  The  measurement  vector,  measurement  covariance  and  particle  filter  state  vector  used  for  all 
sensor  types  are  defined  as 

Ik  =izr,k’zQ,k]' ,  Rk=diag{rj] \  a\ '  ),  xk  =[rk,Qk,kk,Qk]' .  (4.2.1) 

^ r,k  * 6,k 

Using  (4.2.1),  a  generalized  likelihood  function  can  be  defined  as  the  product  of  the  individual  range/bearing 
likelihoods  for  a  given  sensor  type  as  follows: 

1.  Range-Only  Sensor  -  p(zk  \  xh )  =  p(z  k  \  xk)p(zQk  \  xk)  =  N(z  k,al  )U (6 min,9 max) 

^ r,k 

2.  Bearing-Only  Sensor-  p(zk  \  xk)  =  p(zrk \xk)p(zQk  \  xk)  =  U(rmin,rmax)N(zQk,a2_^)  (4.2.2) 


(4.1.1) 

(4.1.2) 
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3.  Range/Bearing  Sensor  -  p(zk  \xk)  =  p(z  k  \  xk)p(zek  \xk)  =  N(z  k,a2  )N(z0k,a2  ) 

^ r,k  ^Q./c 

2  ....  ? 
where  N{z,a“)  corresponds  to  a  normal  distribution  with  mean  z  and  variance  a,  and  U  (amin,amwc ) 

corresponds  to  a  uniform  distribution  over  the  region  amin  to  amax. 


The  generalized  transition  density  used  for  this  work  is  based  on  a  zero-mean,  first-order  Markov  model  [2 1  ] 
and  is  defined  by  the  following  state  transition  matrix  and  process  noise  covariance 


1  0  zr 

(]-e~At/xr 

) 

0 

0  1 

0 

Te  (1- 

.g-AGx0) 

0  0 

g-A  t/zr 

0 

’ 

0  0 

0 

e 

-At  /  Tg 

2a  2  /  zr 

0 

0 

0 

A t3  /3 

0 

At1  /  2 

0 

0 

2gq  /  re 

0 

0 

0 

At3  /3 

0 

At 2  /  2 

0 

0 

2a].  /  zr 

0 

A?2  /  2 

0 

At 

0 

0 

0 

0 

2g@  /  zq_ 

0 

At2  / 2 

0 

At 

where  zr,  zg,  ay,  and  og  are  tunable  parameters  related  to  the  expected  target  motion  in  range/bearing  space  - 
each  of  which  can  assume  different  values  depending  upon  the  amount  and  type  of  fusion  which  has  occurred 
for  a  specific  posterior. 

Applying  the  generalized  likelihood  function  and  transition  density  defined  in  (4.2.2)  and  (4.2.3)  to  (4.1.1) 
and  (4.1.2),  the  particle  filter  is  observed  to  provide  a  consistent  method  of  local  posterior  estimation  for  the 
multi-source  sensor  types  evaluated  in  this  work. 


4.3  Parameterizing  Local  Posteriors 

Communication  of  information  is  paramount  for  any  form  of  distributed  data  fusion.  Consequently,  a 
fundamental  requirement  for  decentralized  fusion  is  efficient  scaling  of  communicated  information  relative  to 
any  growth  in  the  number  of  targets  or  region  of  coverage.  To  date,  the  majority  of  methods  investigated  for 
communicating  information  in  decentralized  architectures  have  been  focused  on  particle  set  optimization 
[7,8],  Gaussian  mixture  models  or  Parzen  density  estimators  [9,10]  or  some  combination  of  the  two 
[1 1,12,13].  All  of  these  methods,  however,  are  prone  to  scalability  concerns  as  the  number  of  particles,  GMM 
components  or  Parzen  components  required  to  adequately  estimate  and/or  reconstruct  posteriors  of  interest 
necessitate  considerable  increase  as  the  number  of  targets  or  region  of  coverage  escalate.  Our  solution  to  this 
scalability  concern  is  rooted  in  high-order  statistics  and  the  generalized  lambda  distribution  (GLD)  -  a  four- 
parameter  probability  distribution  whose  density  function  can  assume  a  wide  variety  of  uni-modal,  symmetric 
and/or  asymmetric  shapes  [22, 23]3. 


3  The  uni-modal  limitation  of  the  generalized  lambda  distribution  prevents  true,  decentralized  multi-modal  estimation.  However,  for 
applications  such  as  target  tracking  where  a  uni-modal  assumption  for  the  marginal  distributions  of  interest  may  not  be  too 
erroneous,  the  generalized  lambda  distribution  provides  an  intermediate  solution  to  the  decentralized  scalability  concern  until  a 
comparable  multi-modal  version  can  be  developed. 
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Derived  from  Tukey’s  original  lambda  distribution  defined  by  the  quantile  function 


Q(u)  =  - 


uL  -  (1  -  u ) 
log(u ) 


X 


(l-u) 


,  A,  ^  0 
,L  =  0 


(4.3.1) 


for  0  <  u  <  1  [24],  the  generalized  lambda  distribution  parameterizes  (4.3.1)  as  follows 


Q(u )  —  H — 


-(1  -u)X4 
^2 


(4.3.2) 


where  A]  acts  as  a  location  parameter,  acts  as  a  scale  parameter,  and  the  combination  of  A,3  and  a4  jointly 
capture  the  shape  of  the  empirical  distribution. 


In  order  to  describe  a  probability  distribution  using  the  generalized  lambda  distribution,  the  ^-parameters  are 
typically  matched  to  the  first  four  moments  of  a  desired  data  distribution  using  a  methodology  similar  to  the 
one  outlined  in  [22].  Figure  6  illustrates  the  ability  of  the  generalized  lambda  distribution  to  accurately 
estimate  various  uni-modal  data  distributions.  The  first  four  moments  of  the  original  data  distribution  and 
those  of  the  GLD  distribution  are  annotated  on  each  plot.  From  Figure  6,  the  mean  percent  moment  error  is 
observed  to  be  less  than  1%  for  all  GLD  estimated  distributions. 


Normal  (100,20) 


Uniform  (100,200) 


Beta  (10,1)  Gamma  (5,10) 

Figure  6:  GLD  Estimation  of  Various  Uni-Modal  Distributions 
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Table  1  -  taken  from  [11,12]  and  slightly  modified  -  illustrates  the  bandwidth  requirements  for 
communicating: 

•  a  particle  set 

•  a  GMM  approximation  of  a  particle  set 

•  a  Parzen  density  approximation  of  a  particle  set 

•  a  GLD  approximation  of  a  particle  set 

In  [11,12]  it  was  stated  that  20  GMM  components  and  50  Parzen  components  were  required  to  adequately 
approximate  a  set  of  2000  particles.  Extending  these  parameters  as  hard  limits  for  accurate  posterior 
representation  and  noting  that  the  GLD  requires  4  components  per  dimension  of  the  state  vector,  Table  1 
illustrates  the  number  of  posteriors  that  can  be  accurately  estimated  for  a  given  state  vector  dimensionality 
and  a  given  amount  of  available  bandwidth  using  the  following  formula: 

„  T  _  .  ,  Available  Bandwidth  (floats) 

Number  ot  Posteriors  Accurately  Communicated  =  - - — - -  . 

Required  Bandwidth  (floats)  Y  Required  Components  | 

Component  J!  Accurate  Posterior  Representation  J 

From  Table  1  it  can  be  seen  that  the  generalized  lambda  distribution  provides  a  significantly  more  compact 
method  of  posterior  representation  as  compared  to  all  the  other  methods  considered. 


Table  1:  Bandwidth  Requirements 


Representation 

Dimension 

Available 

Bandwidth 

(floats) 

Required 
Bandwidth 
(floats)  per 
Component 

Required  Components 
per  Accurate  Posterior 
Representation 

Number  of 
Posteriors 
Accurately 
Communicated 

Particle 

4 

500 

1 

2000 

0 

GMM 

4 

500 

14 

20 

1 

Parzen 

4 

500 

5 

50 

2 

GLD 

4 

500 

1 

16 

31 

Particle 

6 

500 

1 

2000 

0 

GMM 

6 

500 

27 

20 

0 

Parzen 

6 

500 

7 

50 

1 

GLD 

6 

500 

1 

24 

20 

4.4  Reconstructing  Local  Posteriors 

One  advantage  of  using  a  quantile  function  like  (4.3.2)  to  describe  a  distribution  lies  in  its  ability  to  efficiently 
generate  random  variates.  For  instance,  if  Q  is  the  quantile  function  for  a  specific  probability  distribution, 
then  through  the  use  of  the  inverse  transformation  method, 

X  =  Q{u),  u  e  (0,1),  (4.4.1) 
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random  variates,  X,  with  quantile  function,  Q,  can  be  generated  with  little  effort.  Consequently,  since  sources 
of  uniform  random  variates  over  the  interval  (0,1)  are  readily  available,  (4.4.1)  provides  an  efficient  method 
for  generating  random  variates  from  distributions  whose  quantile  functions  are  known  and  computationally 
tractable.  As  a  result,  and  as  demonstrated  in  Figure  6,  this  property  coupled  with  (4.3.2)  allows  for  efficient 
and  accurate  reconstruction  of  local  posteriors  communicated  throughout  a  sensor  network  via  a  set  of  GLD 
parameters. 


4.5  Estimating  Global  Posteriors 

After  the  local  posteriors  have  been  parameterized,  communicated  and  reconstructed  at  any  receiving 
sensor(s),  the  number  of  targets  and  corresponding  positions  are  estimated  using  an  iterative,  probabilistic 
approach.  All  posteriors  at  a  given  sensor  are  then  associated  with  the  estimated  target  information,  and 
corresponding  global  posteriors  are  determined  by  fusing  the  associated  local  posteriors  via  multiplication  of 
their  joint  probability  densities. 

Joint  probability  densities  are  typically  estimated  using  histogram  methods,  however,  such  methods  often 
introduce  undesirable  “artifact  patterns”  which  hamper  the  global  fusion  process  by  introducing  local  extrema 
and  obscuring  the  true  global  optimum  [25].  Furthermore,  these  artifacts  can  be  amplified  when  the  data  is 
sparsely  distributed  as  is  often  the  case  for  particle  filter  applications.  In  an  attempt  to  avoid  such  artifacts, 
interpolation  based  methods  are  commonly  used,  however,  such  methods  do  not  always  remove  the  initial 
artifacts  of  the  histogram  and  frequently  introduce  new  ones.  Our  solution  to  minimizing  histogram-based 
artifacts  is  based  on  an  approximation  technique  referred  to  as  generalized  partial  volume  estimation  (GPVE). 

GPVE  was  originally  developed  for  medical  and  remote  sensing  image  registration  applications  [25,26,27]. 
The  problem  of  estimating  the  joint  probability  density  of  two  uni-dimensional  data  sequences,  however, 
provides  a  slight  alteration  to  the  image  registration  problem.  Consequently,  in  order  to  apply  the  basic  GPVE 
principles  to  the  application  of  estimating  joint  histograms  from  1-D  data  sequences,  the  original  GPVE 
histogram  estimation  procedure  outlined  in  [25,26,27]  needs  to  be  reformulated  using  the  following  mapping: 


X*  :  x  —>  X*(x),  xg  X 
Y*  :  y  — »  7*  (y),  y  gY 


(4.5.1) 


where  X  is  the  discrete  domain  of  the  data  sequence  X*  and  7  is  the  discrete  domain  of  the  data  sequence  7*. 
In  addition,  the  values  X  (x)  and  7  (y)  represent  the  actual  data  values  of  the  two  data  sequences  at  the 
corresponding  indices  x=x„  y=y/  in  terms  of  the  sample  spacing.  Defining  a  real-valued  kernel  satisfying  the 
constraints 


1.  /(«)>  0,  and 

9  Vtn  i  (4'5‘2) 

n=- oo 

where  nsZ  and  Z  is  the  set  of  all  integers,  for  each  x, e V,  y7e7,  the  joint  histogram  of  X*  and  7*  can  be 
updated  in  the  following  manner: 

h(X*(xj  +  p),Y*  (y  j  +q))<^h(X*(xi  +  p),Y*  (y  j  +  q))  +  fx*  (p)  '  fr*  (<?)  Yp,?  e  Z.  (4.5.3) 
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From  (4.5.3)  it  can  be  seen  that  the  increment  of  the  joint  histogram  is  represented  in  terms  of  the  kernel 
functions  corresponding  to  each  data  sequence.  Furthermore,  the  first  constraint  in  (4.5.2)  ensures  that  the 
increments  are  non-negative  while  the  second  constraint  in  (4.5.2)  ensures  that  the  sum  of  the  updated 
amounts  are  normalized  to  one  for  each  corresponding  pair  of  indices  x,  in  X*  and  y,  in  Y*. 


The  appropriate  selection  of  a  suitable  kernel  for  use  in  (4.5.3)  can  depend  on  many  factors,  some  of  which 
may  be  the  histogram  grid  parameters  or  the  sparseness  of  the  data  in  X  relative  to  Y  .  Figure  7  illustrates 
various  kernels  of  the  form / ~  V(0,g2),  all  of  which  satisfy  (4.5.2). 


N(0,15)  N(0,5)  N(0,1)  N(0,0.05) 


n  n  n  n 


Figure  7:  Various  GPVE  Kernels 


Figure  8  illustrates  a  range-only  particle  distribution  with  standard,  cubic-interpolated  and  GPVE  histogram 
estimation  -  all  independently  normalized.  From  this  figure  it  can  be  seen  that  the  standard  histogram 
generates  a  significant  amount  of  distortion  in  the  estimated  distribution  due  to  the  sparseness  of  the  particle 
data.  Consequently,  a  poor  representation  for  the  range-only  particle  distribution  results.  The  interpolated 
histogram  minimizes  some  of  the  distortion  observed  in  the  standard  approach,  however,  a  substantial  amount 
of  distortion  still  remains.  Finally,  the  GPVE  histogram  is  observed  to  accurately  reflect  the  true  distribution 
of  the  range-only  particle  distribution. 


Particle  Distribution 


Histogram 


Northing  (meters)  -100  -100  Easting  (meters) 


Interpolated  Histogram 


Northing  (meters)  -1 00  -1 00  Easting  (meters) 


Generalized  Partial  Volume  Estimation 
x  10'3 


Figure  8:  Joint  Histogram  Estimation  of  Range-Only  Particle  Distribution 


RTO-MP-SET-107 


13  - 13 


UNCLASSIFIED/UNLIMITED 


UNCLASSIFIED/UNLIMITED 


Efficient  Multi-Source  Data  Fusion  for  Decentralized  Sensor  Networks 


ORGANIZATION 


5.0  PROBABILISTIC  DDF  -  SIMULATED  PERFORMANCE  RESULTS 

Figure  9  illustrates  the  target  path  and  sensor  configuration  for  a  simulated  scenario  consisting  of  two  targets, 
three  range-only  sensors  and  one  bearing-only  sensor.  In  this  simulation,  each  target  traverses  around  the 
target  path  counter-clockwise  one  complete  revolution.  All  local  posteriors  are  estimated  using  2000 
particles.  The  sensor  likelihood  parameters  used  in  (4.2.2)  are  given  as  follows: 

1 .  Range-Only  Sensors 

<7Z  =50  meters  downrange  error  at  500  meters 

6  min  =  0  degrees 
6  max  =360  degrees 

2.  Bearing-Only  Sensor 

aZg  -  5  degrees  (90  meters  crossrange  error  at  500  meters) 
rmm  =0  meters 
'max  =  2000  meters 


-400  -200  0  200  400  600  800 

Easting  (meters) 


Figure  9:  Simulation  Scenario 


Figure  10  illustrates  snapshots  of  local  posterior  (top  row)  and  global  posterior  (bottom  row)  particle 
distributions  -  color-coded  relative  to  a  specific  sensor  -  for  no  fusion,  3  fusions,  6  fusions  and  10  fusions. 
From  this  figure  it  can  be  seen  that  the  fused  global  posteriors  for  all  sensors  quickly  collapse  to  the  associated 
target  of  interest. 
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Figure  10:  Particle  Distribution  Snapshots 


Figure  1 1  illustrates  the  norm  of  the  position  error  between  each  target  and  the  mean  of  the  corresponding 
associated  particle  set  at  the  bearing-only  sensor4.  From  Figure  11  it  is  interesting  to  observe  that  certain 
regions  exist  where  the  sensor  modalities  are  not  as  complementary  as  other  regions  (i.e.,  approximately  60, 
130,  300  and  390  seconds).  These  regions  of  concern  can  be  compensated  for,  however,  by  either  adding 
more  sensors,  reconfiguring  the  sensor  modalities  or  rearranging  the  sensor  placement.  Regardless,  from 
Figures  10  and  11  it  can  be  seen  that  the  probabilistic  architecture  outlined  in  this  work  is  capable  of 
accurately  and  efficiently  fusing  multi-source  sensor  data  in  non-linear/non-Gaussian  environments. 
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Figure  11:  Position  Error 


4  The  error  statistics  for  the  range-only  sensors  are  nearly  identical  to  the  error  statistics  of  the  bearing-only  sensor  illustrated  in 
Figure  11. 
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6.0  SUMMARY  AND  FUTURE  WORK 

A  highly-scalable,  Bayesian  approach  to  the  problem  of  performing  multi-source  data  fusion  and  target 
tracking  in  decentralized  sensor  networks  has  been  presented.  The  proposed  probabilistic  approach  provides 
full  decentralization;  seamless  integration  and  efficient  fusion  of  multi-source  sensor  data,  and  an 
unprecedented  scaling  capability.  Two  primary  limitations  of  the  proposed  architecture,  however,  lie  in  its 
current  inability  to  parameterize  multi-modal  distributions  for  information  communication  and  its  lack  of  a 
method  for  handling  common  information  during  fusion.  In  order  to  maintain  efficient  scalability,  it  is 
believed  that  methods  of  density  parameterization  and  fusion  should  continue  to  be  rooted  in  techniques 
utilizing  high-order  statistics.  Consequently,  a  possible  solution  for  parameterizing  multi-modal  distributions 
via  high-order  statistics  could  be  found  through  quantile  or  copula  mixtures.  In  addition,  a  possible  solution 
for  removing  common  information  during  fusion  could  also  be  found  through  a  quantile  or  copula  based 
method  approaching  a  similar  technique  as  outlined  in  [12]. 
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Decentralized  Data  Fusion  -  What  Is  It? 


-  Efficient  method  of  data  fusion  for  distributed  sensor  networks. 


-  Each  node  in  distributed  network  is  a  fusion  node. 


-  Each  fusion  nodes  possesses: 

-  Processing  capability  for 
generating  local  image  of 
fusion  environment. 

-  Capacity  to  share  local 
information  with  neighboring 
nodes  without  requiring 
knowledge  of  network 
topology. 
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Decentralized  Data  Fusion  -  Benefits 


-  Robust: 

-  No  single  point  of  failure. 

-  Modular: 

-  Nodes  are  able  to  dynamically  enter  and/or  leave  network. 

-  Scalable: 

-  Scales  directly  to  number  of  targets  (independent  to  number  of  nodes 
and  measurements). 

-  Channel  management  governs  information  flow  providing  bandwidth 
conservation. 

-  “Common  Map”: 

-  Autonomous  (non-GPS)  localization  and  mapping  ensures  coherent 
picture  of  fusion  environment  across  network  regardless  of  local 
observability  limitations. 
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Linear  Decentralized  Data  Fusion 

-  Based  on  formal  information  metrics  for  quantifying,  communicating 
and  assimilating  sensory  data  in  distributed  sensor  networks. 

-  Given  standard,  linear  discrete  model  of  fusion  environment 

xk  =  Fkxk- 1  +  Gkwk 
zk  =Hkxk  +vk’ 

local  estimates  of  fusion  environment  maintained  at  each  node  using 
Kalman  or  Information  filter: 


Kalman  Filter  Form 

x(k\k),P(k\k) 


Information  Filter  Form 

y(k\k),Y(k\k) 


Kalman/Information  Relationship 


y(k\k)  =  P(k\k)A  x(k\k) 
Y{k\k)  =  P{k\kf 
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-  Information  form  of  Kalman  filter  provides  simple  additive  framework 
for  fusing  target  tracks 

y(k\k)+  =ya{k\k)+yb(k\k)-yab(k\k-\) 

Y(k  |  k)+  =Ya(k\  k)+Yb(k  \  k)-Yab(k  \k-\) 

where  and  Yab  correspond  to  the  “common  information”  between 
nodes  a  and  b. 

-  If  not  handled  correctly,  re-fusion  of  common-information  ultimately 
degrades  fusion  quality. 
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Channel  Filters 


-  Probability  distribution  conditioned  on  common  information  held 
between  a  pair  of  nodes  -  p(xk  \  zka  nZkb ). 

-  Allows  global  posteriors  to  be  estimated  through  communication  of 
local  posteriors  instead  of  likelihood  or  measurement  information. 

k  _  p(xk\Zk)p(xk\Zb )  y(k  |  kf  =ya{k\k)  +  yb(k  \  k)  -  yab(k  \k-\) 
k  P(xk\ZkanZkb )  Y{k\kf  =  Ya(k\k)  +  Yb(k\k)-Yab(k\k-\) 

-  Minimizes  bandwidth  -  fused  estimates  scale  better  than  likelihood  or 
measurement  information. 

-  Constrained  to  unique,  single-connected  (i.e.,  “tree”)  networks. 

-  For  arbitrary  networks,  Covariance  Intersection  (Cl)  algorithm  can  be 
used  to  estimate  common  information  between  nodes. 
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-  Provides  solution  to  problem  of  combining  two  random  vectors  when 
correlation  between  vectors  is  unknown. 

-  Convex  combination  of  information  vectors  and  matrices 

y(k  |  k)+  =  coya(k  \  k)  +  (1  -  co)yb(k  \  k) 

Y(k  |  k)+  =coYa(k\k)  +  (\- co)Yb(k  \  k) 


- Covariance  1 

-  - Covariance  2 

- Covariance  Update  (Cl) 

4 

3 

2 

1 

- Covariance  1 

- - Covariance  2 

- Covariance  Update  (Optimal) 
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\ 

V  ; 

a) 

£  0 
cc 

1 

=?  -1 

-2 

-  Y^~ 

- 

-3 

- 

-4 

x-variance  (m2)  x-variance  (m2) 


Covariance  intersection  guarantees  consistent  estimates  at  expense  of  optimality. 
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Probabilistic  Decentralized  Data  Fusion 


-  Linear  DDF  effectively  used  for  number  of  tracking,  classification  and 
Simultaneous  Localization  and  Mapping  (SLAM)  applications. 

-  Existing  linear  DDF  implementations  restricted  to  uni-modal/uni¬ 
source  sensor  networks  using  Gaussian  based  approaches. 

-  Growing  interest  to  employ  multi-modal/multi-source  sensing 
capabilities  for  improved  “sensor  orthogonality”. 

-  Multi-modal  sensors  potentially  exhibit  observation  and/or  process 
non-linearities  along  with  non-Gaussian  distributions. 

-  Generalized  and  scalable  method  of  decentralized  data  fusion 
required  to  properly  handle  sophisticated  multi-modal  sensor  types. 
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Probabilistic  Decentralized  Data  Fusion  Architecture 


-  Generalized,  Bayesian  approach  to  performing  multi-source  data 
fusion  and  target  tracking  in  decentralized  sensor  networks: 
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-  Point  mass  representations  of  probability. 

-  Approximate  optimal  Bayesian  recursion  without  requiring  model  or 
noise  distribution  constraints. 

-  Posterior  density  of  interest  represented  by  set  of  weighted  particles. 

-  Overall  estimate  of  true  posterior  density  obtained  from  particle 
ensemble. 


Multi-Source  Sensor  Particle  Distributions  Multi-Source  Sensor  Probability  Distributions 
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Estimating  Local  Posteriors  -  Generalized  Likelihood  Function 


-  Generalized  likelihood  function  defined  as  product  of  individual 
range/bearing  likelihoods  for  specific  sensor: 

-  Range-only: 

p(z_fc  |  X/{ )  =  p(zr,k  I  =£ k  )p(zQ,k  I  ^k )  =  ^(zr,k  > az  )U{@min> @max ) 

r,k 

-  Bearing-only: 

=  U*b(ze,/t  UA:)  =  c/(o«i«<w)Ar(ze^-0?  ) 

z0,£ 

-  Range/Bearing: 

PUk \xk)  =  P(zr,k \xk)p(zQ,k  \xk)=N{zrk,a2  )N(zQk,o2  ) 

Zr,k  ZQ,k 

where 

N{z,o2)  -  normal  distribution  with  mean  z  and  variance  a2. 
U(amin,am J  =  uniform  distribution  over  region  amin  to  amax. 
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Estimating  Local  Posteriors  -  Generalized  Transition  Density 

-  Generalized  transition  density  based  on  zero-mean,  first-order 
Markov  model: 


a>  = 


1  0  t  {l-e 


-At  /  T, 


0  1 
0  0 
0  0 


2  a  2  /  Tr 


Qk  = 
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Given  the  appropriate  likelihood  functions  and  transition  densities , 
the  particle  filter  provides  a  consistent  method  of  local  posterior 
estimation  for  the  multi-source  sensor  types  evaluated  in  this  work. 
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Communicating  Information  In  Distributed  Networks 


-  Communication  of  local  information  is  paramount  for  any  form  of 
distributed  data  fusion. 

-  Efficient  scaling  of  communicated  information  relative  to  growth  in 
number  of  targets  or  region  of  coverage  is  fundamental. 

-  Typical  methods  currently  used  to  communicate  information  in 
probabilistic  decentralized  architectures: 

1 .  Particle  set  optimization. 

2.  Gaussian  mixture  models  (GMMs). 

3.  Parzen  density  estimators. 

-  Such  methods  are  prone  to  scalability  concerns. 

-  Our  solution  to  this  scalability  concern  is  rooted  in  high-order 
statistics  and  the  generalized  lambda  distribution  (GLD). 
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Generalized  Lambda  Distribution 


-  Four-parameter  probability  distribution  whose  density  function 
assumes  a  wide  variety  of  uni-modal,  symmetric  and/or  asymmetric 
shapes. 

-  Parameterizes  a  data  distribution  by  matching  first  four  moments  of 
distribution  to  ^-parameters  of  following  quantile 


Q(u)~  +  — 


-(1-k)X4 

h 


where  acts  as  a  location  parameter,  X2  acts  as  a  scale  parameter 
and  A.3,A,4  jointly  capture  shape  of  empirical  distribution. 


13-15 


BAE  SYSTEMS 


Generalized  Lambda  Distribution  (continued) 

-  Illustrations  of  GLD  estimating  various  uni-modal  data  distributions. 

-  Mean  percent  moment  error  <  1%  for  all  GLD  estimated  distributions. 


Normal  (100,20) 


Uniform  (100,200) 


Beta  (10,1) 


Gamma  (5,10) 


The  generalized  lambda  distribution  is  capable  of  accurately  and 
efficiently  estimating  a  wide-variety  of  uni-modal  data  distributions. 
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-  Comparison  of  bandwidth  requirements  for  various  methods  of 
information  communication: 

-  Particle  Sets: 

Assumption  -  2000  particles  required  for  accurate  posterior  representation. 

-  Gaussian  Mixture  Models: 

Assumption  -  20  GMM  components  required  for  accurate  posterior  representation. 

-  Parzen  Density  Estimators: 

Assumption  -  50  Parzen  components  required  for  accurate  posterior  representation. 

-  Generalized  Lambda  Distribution: 

4  components  per  state  dimension  required  for  accurate  posterior  representation. 
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Communication  Bandwidth  Comparison  (continued) 


-  Number  of  posteriors  accurately  communicated  using  various  methods 
of  information  communication  with  corresponding  assumptions: 


Representation 

State 

Dimension 

Available 

Bandwidth 

(floats) 

Required  Bandwidth 
(floats) 

per  Component 

Required  Components 
per  Accurate 
Posterior  Representation 

Number  of  Posteriors 
Accurately 
Communicated 

Particle 

4 

500 

1 

2000 

0 

GMM 

4 

500 

14 

20 

1 

Parzen 

4 

500 

5 

50 

2 

GLD 

4 

500 

1 

16 

31 

Particle 

6 

500 

1 

2000 

0 

GMM 

6 

500 

27 

20 

0 

Parzen 

6 

500 

7 

50 

1 

GLD 

6 

500 

1 

24 

20 

Number  of  Posteriors  Accurately  Communicated  = 


Available  Bandwidth  (floats) 


Required  Bandwidth  (floats) 

Required  Components 

v  Component  , 

Accurate  Posterior  Representation , 

The  generalized  lambda  distribution  provides  a  significantly  more 
compact  method  of  posterior  representation  over  all  other  methods. 
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Reconstructing  Local  Posteriors 


-  Given  a  quantile  function  for  a  specific  probability  distribution, 
random  variates,  X,  with  quantile,  Q,  can  be  generated  with  little 
effort  using  the  inverse  transformation  method 

X-Q(u),  ug(0,1). 


-  Sources  of  uniform  random  variates  over  the  interval  (0,1)  are 
readily  available. 

-  Inverse  transformation  provides  efficient  and  accurate  method  for 
reconstructing  local  posteriors  communicated  throughout  sensor 
network  via  GLD  parameters. 
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-  After  local  posteriors  have  been  received  and  reconstructed  at  a 
given  sensor,  the  number  of  targets  and  corresponding  positions  are 
estimated  using  an  iterative,  probabilistic  approach. 

-  Local  posteriors  are  associated  with  approximated  target  information. 


-  Global  posteriors  are  estimated  by  multiplying  joint  probability 
densities  of  associated  local  posteriors. 


•1000  -800  -600  -400  -200  0  200  400  BOO  800  1000  -1000  -800  -600  -400  -200  0  200  400  600  800  1000  -1000  -800  -600  -400  -200  0  200  400  600  800  1000 

Easting  (meters)  Easting  (meters)  Easting  (meters) 
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Simulated  Performance  Results 


-  Target  path  and  sensor  configuration  for  simulated  scenario  consisting 
of  two  targets,  three  range-only  sensors  and  one  bearing-only  sensor. 

Road  O  Target  Range-Only  Sensor  Bearing-Only  Sensor  Comm  Link 
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Simulated  Performance  Results  -  Parameters 

-  Local  posteriors  estimated  using  2000  particles. 

-  Range-only  sensor  likelihood  parameters: 

az  =50  meters  downrange  error  at  500  meters 

0mm  =  0  degrees 
Omax  =  360  degrees 

-  Bearing-only  sensor  likelihood  parameters: 

crZ0  =  5  degrees  (90  meters  crossrange  error  at  500  meters) 

rmin  =  0  meters 
rmax  =  2000  meters 
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Simulated  Performance  Results  -  Posterior  Illustrations 


-  Snapshots  of  local  and  global  posterior  particle  distributions. 
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Simulated  Performance  Results  -  Position  Error 

-  Norm  of  position  error  between  each  target  and  mean  of  corresponding 
associated  particle  set. 
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Decentralized  probabilistic  architecture  is  capable  of  accurately  and  efficiently 
fusing  multi-source  sensor  data  in  non-linear/non-Gaussian  environments. 
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Summary 


-  A  highly-scalable,  Bayesian  approach  to  multi-source  data  fusion 
and  target  tracking  in  decentralized  sensor  networks  has  been 
presented. 

-  Proposed  probabilistic  approach  provides: 

1.  Full  decentralization. 

2.  Seamless  integration  and  efficient  fusion  of  multi-source  sensor  data. 

3.  Unprecedented  scaling  capability. 

-  Simulated  results  illustrate  capability  of  proposed  probabilistic 
architecture  to  accurately  and  efficiently  fuse  multi-source  sensor 
data  in  complex,  non-linear/non-Gaussian  environments. 

-  Current  Research: 

1.  Extending  GLD  uni-modal  parameterization  to  multi-modal  capability. 

2.  Developing  methods  to  estimate  and  remove  common  information 
during  fusion. 
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